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We study a model of quasiparticles on a two-dimensional square lattice coupled to Gaussian dis- 
tributed dynamical molecular fields. We consider two types of such fields, a vector molecular field 
that couples to the quasiparticle spin-density and a scalar field coupled to the quasiparticle number 
density. The model describes quasiparticles coupled to spin or charge fluctuations, and is solved by 
a Monte Carlo sampling of the molecular field distributions. The nonperturbative solution is com- 
pared to various approximations based on diagrammatic perturbation theory. When the molecular 
field correlations are sufficiently weak, the diagrammatic calculations capture the qualitative aspects 
of the quasiparticle spectrum. For a range of model parameters near the magnetic boundary, we 
find that the quasiparticle spectrum is qualitatively different from that of a Fermi liquid, in that 
it shows a double peak structure, and that the diagrammatic approximations we consider fail to 
reproduce, even qualitatively, the nonperturbative results of the Monte Carlo calculations. This sug- 
gests that the magnetic pseudogap induced by a coupling to antiferromagnetic spin-fluctuations and 
the spin-splitting of the quasiparticle peak induced by a coupling to ferromagnetic spin-fluctuations 
lie beyond diagrammatic perturbation theory. While a pseudogap opens when quasiparticles are 
coupled to antiferromagnetic fluctuations, such a pseudogap is not observed in the corresponding 
charge-fluctuation case for the range of parameters studied, where vertex corrections are found to 
effectively reduce the strength of the interaction. This suggests that one has to be closer to the 
border of long-range order to observe pseudogap effects in the charge-fluctuation case than for a 
spin-fluctuation induced interaction under otherwise similar conditions. The diagrammatic approx- 
imations that contain flrst order vertex corrections show the enhancement of the spin-fluctuation 
induced interaction and the suppression of the effective interaction in the charge-fluctuation case. 
However, for the range of model parameters considered here, the multiple spin or charge-fluctuation 
exchange processes not included in the diagrammatic approximations considered are found to be 
important, especially for quasiparticles coupled to charge fluctuations. 

PACS numbers: PACS Nos. 71.27.-|-a 



I. INTRODUCTION 



The concept of elementary excitations and the di- 
agrammatic perturbation-theoretic methods borrowed 
from quantum field theory have given us, over the past 
decades, many powerful insights into the behavior of 
materials. In a number of cases, however, these con- 
cepts and methods don't seem to work. In previous 
paperai*^., we presented results on a nonperturbative ex- 
tension of the magnetic interaction model, which had 
until then been extensively used in the context of dia- 
grammatic approaches. These latter applications were 
successful in many respects: in the Eliashberg approx- 
imation, the magnetic interaction model correctly an- 
ticipated the pairing symmetry of the Cooper state in 
the copper oxide superconductorsi and is consistent with 
spin-triplet p-wave pairing in superfluid '^He [for a re- 
cent review see, e.g., ref.*]. One also gets the correct 
order of magnitude of the superconducting and super- 
fluid transition temperature Tc when the model param- 
eters are inferred from experiments in the normal state 
of the above systems. However, in Ref.^ it was found 
that when the model was treated nonperturbatively and 
one approached the border of magnetic long-range order, 
the quasiparticle spectrum showed qualitative changes 
not captured by the Eliashberg approximation. In Refi, 



we raised the possibility that these qualitative changes, 
namely the opening of a pseudogap in the quasiparti- 
cle spectrum, were intrinsically nonperturbative in na- 
ture. In this paper, we examine this possibility by com- 
paring the nonperturbative results to various kinds of 
perturbation-theoretic approximations . 

The paper is organized as follows. In the next section 
we describe the model as well as the various perturbation- 
theoretic approximations to be compared to the Monte 
Carlo calculations. Section III contains the results of the 
nonperturbative and diagrammatic calculations. Section 
IV contains a discussion of the results and finally we give 
a summary and outlook 



II. MODEL 

The model and its motivation have been extensively 
discussed in Refi. Here we only give the definitions rel- 
evant to the present discussion. We consider particles on 
a two-dimensional square lattice whose Hamiltonian in 
the absence of interactions is 
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where tij is the tight-binding hopping matrix, fi the 
chemical potential and tpj^, ipia respectively create and 
annihilate a fermion of spin orientation a at site i. We 
take tij = t if sites i and j are nearest neighbors and 
tij — t' if sites i and j are next-nearest neighbors. 

To introduce interactions between the particles, we 
couple them to a dynamical molecular (or Hubbard- 
Stratonovich) field. It is instructive to consider two dif- 
ferent types of molecular fields. In the first instance, we 
consider a vector Hubbard-Stratonovich field that cou- 
ples locally to the fermion spin density. We also con- 
sider the case of a scalar field that couples locally to 
the fermion number density. This case corresponds to a 
coupling to charge-fluctuations or, within the approxima- 
tion we are using here, "Ising"-like magnetic fluctuations 
where only longitudinal modes are present. The Hamil- 
tonians at imaginary time r for particles coupled to the 
fluctuating exchange or scalar dynamical field are then 



hir) = hoir) _ ^^M,(r) . 4jT)a^^i;,^{T) (2) 
h{r) = ho{T)-gJ2'^^ir)4o.ir)4'UT) (3) 



where M,(r) = {M^ {t), Mf (t), Mf (t))'^ and $,(t) 
are the real vector exchange and scalar Hubbard- 
Stratonovich fields respectively, and g the coupling con- 
stant. The reason for the choice of an extra factor l/\/3 
in Eq. Q becomes clear later. 

Since we ignore the self-interactions of the molecular 
fields, their distribution is Gaussian and given bjiS. 



r[M] = iexp ( - 5: i (4) 



Z ^ DM exp 



E 



2a(q, ii^n) 

M(q, ^M(-q,-^i^„)^ 
2a(q, w„) 



(5) 



We consider the case where there is no long-range mag- 
netic or charge order. The average of the dynamical 
molecular fields must then vanish and their Gaussian dis- 
tributions Eqs. H4I6|) are completely determined by their 
variance Q;(q, which we take to be 



5^x(q7«i^n) if M(q, ii/„) or $(q,w„) 



is complex 
if M(q 
is real 



^X(q, -ii^n) if M(q, w„) or $(q, w„) 



where N is the number of allowed wavevectors in the 
Brillouin zone. Then 



M,{q,iiy„)Mj{k,tnr. 



T 



X Sci--kd^^-nji,j (11) 
T 

$(q,iz^„)$(k,il7„)) = — x(q, w„)(5q -k 



(12) 



where < . . . > denotes an average over the probability 
distributions Eq. and Eq. ^ for the vector and scalar 
cases respectively. In order to compare the scalar and 
vector molecular fields, we take the same form for their 
correlation function x(qj ^'^n) parametrize it as in 
Refs.^-*&. In what follows, we set the lattice spacing a to 
unity. For real frequencies, we have 



x(q,w) 



(p — i 



(13) 
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where k and kq are the correlation wavevectors or inverse 
correlation lengths in units of the lattice spacing, with 
and without strong correlations, respectively. Let 



in the case of a vector exchange molecular field and 



Vm = - exp 



^ <I>(q, z;/„)$(-q, -ii^n) j 



Z = / £>$ exp I - ^ 



2a(q, Wn) 

^(q, z;/„)$(-q, -^t^n) 
2a (q, zj/„) 



(7) 



in the case of a scalar Hubbard-Stratonovich field. In 
both cases Vn — 27rnT since the dynamical molecular 
fields are periodic functions in the interval [0,f3 — 1/T]. 
The Fourier transforms of the molecular fields are defined 
as 



Mr{t) = ^M(q,ii/„)exp(^-2[q-R-i/„T]j (8) 
^r{t) = ^$(q,w„)exp('-i[q-R-t/„T]] (9) 



4±2(cos(gj;) +cos{qy)) 



(14) 



We consider commensurate charge fluctuations and an- 
tiferromagnetic spin fluctuations, in which case the pa- 
rameters and ri{q) in Eq. H13I) are defined as 



(15) 
(16) 



where Tq is a characteristic temperature. 

We also consider the case of ferromagnetic spin- 



fiuctuations, where the 
Eq. l(T^ are given by 



parameters q^ and r]{q) in 



q 

viq) 



Taq- 



(17) 
(18) 
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x(q, zf„) is related to the imaginary part of the response 
function /mx(q, a;), Eq. via the spectral represen- 

tation 



X(q,zt'„) = - 



+ 00 



diu /mx(q, uj) 



(19) 



To get x(q, i'^n) to decay as l/v^ as 00, as it should, 

we introduce a cutoff tuo and take /mx(q, oj) = for 

> Wo- A natural choice for the cutoff is luq = ry(q)Ko. 

In our model, the single particle Green's fimction is the 
average over the probability distributions T'iM] (Eq. g))) 
or V[^] (Eq. ©) of the fermion Green's function in a 
dynamical vector or scalar field. 



g{icrT;ja'T') = J DMV[M] G{iaT;ja'T'\[M]){20) 
g{iaT]ja'T') ^ J D^V[<i>]G{iaT;ja'T\[<^]) (21) 



where 



G{taT;ja'T'\[M] or [$]) = -(r,{V^.(T)^],, (r')}) (22) 

is the single particle Green's function in a dynamical 
molecular field and is discussed at length in Ref.^. In 
evaluating expressions Eqs. H2UI21|1 one is summing over 
all Feynman diagrams corresponding to spin or charge- 
fluctuation exchangesiiSil. The diagrammatic expansion 
of the Green's function, Eq. is shown pictorially in 
Fig. 1. Since in our model no virtual fermion loops are 
present, there is no fermion sign problemSi. 

In this paper we compare the results of the Monte 
Carlo simulations to various diagrammatic approxima- 
tions for the same model. We denote by go{p,iuJn) and 
g{p, ioJn) the bare and dressed quasiparticle propagators 
respectively. They are given by 



t/o(p,«w„) = 



1 



(fp - M) 
1 



ILUn - (ep - ^) - i;(p, ILUn) 



(23) 
(24) 



where S(p, iw„) is the quasiparticle self-energy and 
Ep the tight-binding dispersion relation obtained from 
Fourier transforming the hopping matrix in Eq. 
and /i the chemical potential. We consider four approxi- 
mations to the quasiparticle self-energy I](p,zcij„) whose 
diagrammatic representations are shown in Fig. 2. 

In Fig. 2a, the self-energy is approximated by first or- 
der perturbation theory in the exchange of magnetic or 
charge fluctuations and denoted S]^^*(p, iwn). Fig. 2b 
shows the Eliashberg approximation in which the self- 
energy denoted I]^*'^(p, i^n) is given by the first order 
self-consistent (or Brillouin-Wigner) perturbation theory. 



The expressions for S-'^P*(p, iw„) or T,^'"^{p,iuJn) in the 
case of quasiparticles coupled to magnetic or charge fluc- 
tuations are identical (this is the reason for our choice of 
the factor 1/a/3 in Eq. (O) and given by 

X GoiKin,,) (25) 

Ei*'=(p,iu;„) = g^—'^xip-Ki^n-i^ri) 
k,n„ 

X GiKmn) (26) 

where t/o(k, i^ln) and ^(k, zf2„) are the bare and dressed 
quasiparticle Green's functions defined in Eq. (|23|l and 
Eq. (|24|l respectively. Fig. 2c shows the diagrammatic 
expansion corresponding to second order perturbation 
theory and we denote the self-energy corresponding to 
that approximation S^P*(p, iwn). The second order self- 
consistent approximation to the quasiparticle self-energy, 
denoted E^'"^(p, iwn), is shown diagrammatically in Fig. 
2d. The expressions for I]^^'*(p, ia;„) and S^*^(p,iw„) 
now depend on whether the quasiparticles are coupled to 
the vector Hubbard-Stratonovich field (magnetic fluctua- 
tions) or scalar Hubbard-Stratonovich fleld (charge fluc- 
tuations), because vertex corrections in the two cases do 
not have the same coefficient or even the same sign. The 
expressions for S^P*(p,zCi;„) and E^*'^(p,iw„) for quasi- 
particles coupled to magnetic fluctuations are given by 



T,^P\p,iUn) = g^— ^ x(p-k,icj„-if7„) 
X t/o(k,ir2„) 

X go{'k,inn)^^p\p,iLJr,)go{k,in, 
1 



\ / k,0„k',0' 



X Goik, inn)Goik', in'„)xik - k', iflr, - iQ'J 

X go{p-k + k',iuj„-inr, + in'j (27) 

S2^=(p,iw„) = g'^—'^xiP'-Ki^^n-i^n) 

k,0„ 

\ I k,n„k',n^ 

X g(k, za05(k', T-^M^ - k', i^n - 

X g(p-k + k',iw„-il7„ + if7;) (28) 

In the case of the scalar Hubbard-Stratonovich field, or 
coupling to charge fiuctuations, the corresponding ex- 
pressions are 
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\ / k,n„k',n' 



X g{k, in„)g{k',in'jx{'^ ~ k', m,, - ai'j 

X g(p-k + k',it^„-iai + if^^) (30) 



In Eqs. 127I28I29I30|) . go(k,irj„) and g{k,inn) are the 
bare and dressed quasiparticle Green's functions defined 
in Eq. and Eq. (|21|) respectively. 



A. Mass renormalization parameter 

The strength of the couphng to the magnetic or charge 
fiucutations can be parametrized by a dimensionless mass 
renormahzation parameter Xz, which is defined as 



/-rv<^5^^-x(p-p',-) >i^5(p,p' 



< 1 > 



FS{p) 



.(31) 



The Fermi surface averages are given by 



< ■ • ■ >-FS(p) 
< • • • >_FS(p,p') 



X (5(ep - /i)J(ep/ - ^) 



<5(ep - m) (32) 



(33) 



In practice, we compute the Fermi surface average with a 
discrete set of momenta and we replace the delta function 
by a finite temperature expression 



d'^p 
5{ep - n) 



-E 



T 



/p(l-/p) 



(34) 
(35) 



where /p is the Fermi function. Note that y/p(l — /p) 
6{ep-n) as T ^ 0. We have used T = O.lt and N = 128^ 
in all of our calculations. The finite temperature effec- 
tively means that van Hove singularities will be smeared 
out. 

Note that the Fermi surface average that appears in 
Xz, Eq. H31|l plays a role similar to that of a'^F{uj)/uj 
in the case of phonon mediated superconductivity. One 
therefore expects ~ 1 to indicate the crossover be- 
tween weak and strong coupling. 



III. RESULTS 

The quasiparticle dispersion relation for the two- 
dimensional square lattice is obtained from Eq. We 
measure all energies and temperatures in units of the 
nearest-neighbor hopping parameter t. We set the next- 
nearest-neighbor hopping parameter t' = — 0.45t. The 
chemical potential is adjusted so that the electronic band 
filling is n = 0.9. The dimensionless parameters describ- 
ing the molecular field correlations are g^Xo/tj To/t, kq 
and K. We chose a representative value for Kq — 12, and 
set To = 0.67i as in the earlier work^-. For an electronic 
bandwidth of leV, Tq « 1000°K. We only consider one 
value of the coupling constant g'^Xo/t = 2. In the random 
phase approximation, the magnetic instability would be 
obtained for a value of g^Xo/t of the order of 10. We 
consider what happens to the quasiparticle spectrum at 
a fixed temperature T = 0.25t as the inverse correlation 
length K changes, as in Rcf.^. 

All the calculations were done on a 8 by 8 spatial lat- 
tice. In the Monte Carlo calculations we used 41 imag- 
inary time slices, or equivalently 41 Matsubara frequen- 
cies for the molecular fields, M(q, w„) and $(q, w„) 
{i/n = 2TrnT, with n — 0, ±1, . . . , ±20). In the diagram- 
matic calculations, we used between 40 to 60 fermion 
Matsubara frequencies. 

By analytic continuation of the single particle Green's 
function C/(k, r) one can obtain the quasiparticle spec- 
tral function v4(k,a;) = — i/m 5/j(k, w) and the tun- 
neling density of states N{uj) = jj A{k, cu) , where 
Gnlk^Lo) is the retarded single particle Green's function. 
The imaginary time Monte Carlo data is analytically con- 
tinued with the Maximum Entropy method^, using the 
same methodology as in the earlier worki.. We used 10000 
Monte Carlo samples grouped into 100 bins of 100 sam- 
ples each. We always use a flat default model in the 
Maximum Entropy calculations. To provide a fair com- 
parison between diagrammatic and nonperturbative cal- 
culations, one should use the same analytic continuation 
method (with the same parameters) in all cases. There- 
fore, we generated 100 noisy measurements by adding 
Gaussian random noise to the results of the diagram- 
matic calculations and analytically continued G (k, r) = 
Goik, t)+T e-*""^ {G{k, iuJn)-Go{K *a;„)) using the 
Maximum Entropy method as well, with the same de- 
fault model as in the corresponding analytic continua- 
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tion of the Monte Carlo data. The scheme is not per- 
fect, however. While the variance of the Gaussian noise 
added to the diagrammatic Green's functions was cho- 
sen such that the statistical uncertainty of the average 
over the 100 noisy samples was identical to that in the 
corresponding Monte Carlo Green's function, the corre- 
lations in the errors for different values of r present in 
the Monte Carlo results cannot be easily modeled. The 
Gaussian random numbers added to the diagrammatic 
Green's function were therefore taken to be independent 
of each other, and thus the noise in the diagrammatic 
and Monte Carlo Green's functions did not have iden- 
tical statistical properties. In spite of this, the present 
scheme is almost certainly better than the alternatives. 



A. Antiferromagnetic spin-fluctuations 

Figs. 3,4 and 5 show the comparison, for different val- 
ues of K^, between the nonpcrturbative calculations of 
the quasiparticle Green's function fj(k, t), spectral func- 
tion A(k, w) and tunneling density of states N{uj) and 
those obtained from the approximations S^P*(p, icj„), 
S^?'*(p, iujn), S^*'^(p, iujn), and S^*^(p, iujn) to the quasi- 
particle self-energy. 

Fig. 3 shows our results for = 24. For this value 
of the inverse correlation length squared, the mass renor- 
malization parameter Xz ~ 0.05. The coupling to the an- 
tiferromagnetic spin-fluctuations is therefore weak. Not 
surprisingly, the quasiparticle Green's function, spectral 
function and tunneling density of states obtained from 
the various diagrammatic approximation agree well with 
the Monte Carlo results. At k = (tt, 0), the difference 
between the nonpcrturbative Green's function ^((k, r) 
and its diagrammatic approximations is of the order of 
O.OOlt for all values of r. There is virtually no difference 
between the straightforward perturbation-theoretic cal- 
culations of the spectral function and their self-consistent 
counterparts, in first and second order, which is expected 
for weak coupling. Thus the small difference in the spec- 
tral functions A(k, w) at k = (tt, 0), seen in Fig 3b, to the 
extent that they are not an artifact of the analytic contin- 
uation, must come from the vertex corrections. Since the 
first order spectral functions are slightly sharper than the 
second order ones, the first order vertex corrections result 
in an increased spin-fluctuation interaction, as pointed 
out in Refs i^i^i^" . The Monte Carlo spectral function is 
also somewhat broader than the diagrammatic calcula- 
tions, and provided again that it is not an artifact of the 
analytic continuation, this suggests that the higher order 
diagrams lead to a further increase of the spin-fluctuation 
interaction. 

The results for = 4 are shown in Fig. 4. This value 
of gives a mass renormalization parameter Xz ~ 0.45. 
One is now in the intermediate coupling regime. The 
quasiparticle Green's function ^(k, r) and spectral func- 
tion A(k, w) at k = (tt, 0) as well as the tunneling 
density of states obtained from the various diagram- 



matic approximations agree qualitatively with the Monte 
Carlo results. There are, however, noticeable quantita- 
tive differences, not surprisingly much more so than for 

= 24. The largest difference between the Green's 
functions obtained from the diagrammatic approxima- 
tions and the nonpcrturbative calculations is now big- 
ger than the width of the lines and is roughly an or- 
der of magnitude (O.Olt) larger than for = 24, which 
is not unexpected since the mass renormalization pa- 
rameter is also about an order of magnitude greater for 

= 4 than for = 24. With the above caveat regard- 
ing the analytic continuation, one can make a few addi- 
tional remarks. First of all, there is now a difference be- 
tween the straightforward perturbation-theoretic results 
and the self-consistent calculations of the spectral func- 
tion yl(k, w) at k = (tt, 0), both at first and second or- 
der. In particular, the second order self-consistent spec- 
tral function is slightly broader than the first order self- 
consistent one, an indication that the first order vertex 
correction leads to an enhancement of the effective spin- 
fiuctuation interaction, in agreement with Refsi ^i^i'^" . 
The nonpcrturbative A(k = (7r,0),a-') is broader than 
the second order self-consistent result, which would im- 
ply the higher order vertex corrections are further en- 
hancing the magnetic interaction. Note, that the second 
order perturbation-theoretic A(k = {tt,0),uj) is slightly 
broader than its self-consistent counterpart (the dress- 
ing of Green's functions tends to reduce the effect of in- 
teractions) and agrees very well with the Monte Carlo 
result. This may be due to a cancellation of errors (or 
the analytic continuation procedure) since the agreement 
between the nonpcrturbative tunneling density of states 
N{uj) and the second order perturbation-theoretic N{uj) 
is not as good. 

As « 1, the quasiparticle mean free path becomes 
of the order of the magnetic correlation length for some 
wavevectors near the Fermi surface, the quasiparticles 
then can't tell there is no long-range order, and this 
marks the onset of pseudogap behavior^. For ~ 1, 
the mass renormalization parameter is Xz ~ 1.1. One is 
therefore in the strong coupling regime. The results of 
our calculations for = 1 are shown in Fig. 5. The 
developing pseudogap in the spectral function A(k, w) 
at k = (tt, 0) (Fig. 5b) and in the tunnehng density 
of states N{lu) (Fig. 5e) found in the nonpcrturbative 
Monte Carlo calculations is not seen in any of the dia- 
grammatic approximations considered here, which there- 
fore fail qualitatively. Given that one is in the strong 
coupling regime > 1, the breakdown of perturbation 
theory should not come as a surprise. The maximum 
difference in the quasiparticle Green's function tj(k, r) 
between the nonpcrturbative and diagrammatic calcula- 
tions is now of the order O.lt, and hence an order of mag- 
nitude larger than for = 4 and a couple of orders of 
magnitude larger than in the weak coupling regime with 
= 24. It is therefore not suprprising that the quasi- 
particle spectra that give rise to these rather different 
imaginary time Green's functions turn out to show qual- 
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itative differences. Note that for = 1, there is nearly as 
much difference between the perturbation-theoretic and 
self-consistent approximations of the same order as there 
are between calculations of the same type at first and 
second order. 



B. Ferromagnetic spin- fluctuations 

Figs. 6-9 show our results for the quasiparticle Green's 
function C/(k, t), spectral function A(k, tj) and tunneling 
density of states N{ui) for several values of k^. 

We start with = 24, for which the mass renormaliza- 
tion parameter Xz ~ 0.05 for coupling to ferromagnetic 
spin-flucutations. In this weak coupling regime, Fig. 6 
shows that the results of the various diagrammatic cal- 
culations are in good agreement with the Monte Carlo 
results. As in the corresponding antiferromagnetic case, 
at k = (tt, 0), the difference between the nonperturbative 
Green's function t/((k, t) and its diagrammatic approxi- 
mations is of the order of O.OOlt for all values of r. More- 
over, there is virtually no difference between the straight- 
forward perturbation-theoretic calculations of the spec- 
tral function and their self-consistent counterparts, in 
first and second order. Thus the small difference in the 
spectral functions A(k, at k = (7r,0), seen in Fig. 
6b must come from the vertex corrections. Since the 
first order spectral functions are slightly sharper than 
the second order ones, the first order vertex corrections 
result in an increased spin-fluctuation interaction. Given 
the smallness of the difference between the first and sec- 
ond order results and the ill-posed nature of the analytic 
continuation problem, one should take the above remark 
with some degree of caution. In Refi, however, it was 
shown that the increase in the effective interaction in- 
duced by the first order vertex correction is due to the 
spin dependence of the interaction, and thus should oc- 
cur for quasiparticles coupled to either antiferromagnetic 
or ferromagnetic fluctuations. Our analytically contin- 
ued results are at least consistent with this. From Fig. 
6b, one also sees that the Monte Carlo spectral function 
is slightly broader than the first or second order results, 
as in the corresponding antiferromagnetic case. With the 
above caveat on the nature of the analytic continuation 
problem, this would suggest the higher order diagrams 
not included in our perturbation-theoretic approxima- 
tions lead to a further enhancement of the magnetic in- 
teraction, as in the corresponding antiferromagnetic case. 

For k'^ = 4, the mass renormalization parameter Xz ~ 
0.45 and one is therefore in an intermediate coupling 
regime. Fig. 7 shows that for this value of k^, the di- 
agrammatic approximations all qualitatively agree with 
the Monte Carlo results. The quantitative agreement is, 
not surprisingly, not as good as in the weak coupling 
limit with — 24. One notes a number of similari- 
ties between the results of Fig. 6 and the corresponding 
antiferromagnetic case, shown in Fig. 4: (i) the sec- 
ond order perturbation theory results for A(k, lu) give 



the best agreement with the nonperturbative calcula- 
tion, (ii) since the spectral function in either second order 
calculation, which include vertex corrections, is slightly 
broader in uj than the corresponding first order result, 
we conclude that first order vertex corrections lead to 
an enhancement of the effective quasiparticle interaction, 
which is what is expected on the basis of the arguments 
made in Refsi2ii£ (iii) the spectral function obtained by 
Monte Carlo sampling of the Gaussian dynamical molec- 
ular fields is slightly broader than the second order re- 
sults, which to the extent this is not an artifact of the 
Maximum Entropy analytic continuation is an indication 
that higher order spin-fluctuation exchanges not included 
in the diagrammatic approximations considered lead to a 
further enhancement of the effective quasiparticle inter- 
action. 

The dynamical exponent z is larger for ferromagnetic 
than antiferromagnetic spin fluctuations. Hence the ef- 
fective dimension d -I- z = 5 in the ferromagnetic versus 
d + z ~ A in the antiferromagnetic case and the stan- 
dard theory of quantum critical phenomena^ ^ leads one 
to expect weaker corrections for higher effective dimen- 
sions. The perturbative calculations qualitatively fail at 

« 1 in the antiferromagnetic case and on the basis of 
the above arguments one would expect that the break- 
down of perturbation theory in the case of ferromagnetic 
fluctuations, if it happens, would occur for a smaller value 
of or larger values of the mass renormalization param- 
eter Xz- Indeed, at = 1, ~ 1.2 and therefore one 
is in the strong coupling regime. Fig. 8 shows that while 
for this value of the diagrammatic calculations still 
agree qualitatively with the Monte Carlo results, unsur- 
prisingly there are larger quantitative differences than in 
the case — 4 shown in Fig. 7. 

Our results for = 0.25, for which the mass renor- 
malization parameter Xz ~ 2.3 are shown in Fig. 9. 
The spectral function A{k, oj) obtained from the nonper- 
turbative Monte Carlo calculations shows a double peak 
structure. This has been interpreted in Ref.^ as an effec- 
tive spin-splitting of the quasiparticle spectrum induced 
by the local ferromagnetic order. In looking at the evo- 
lution of the spectral function A(k, oj) as is decreased, 
one first sees a broadening of ^(k, ui) and then, the broad 
quasiparticle peak splits into two. The Monte Carlo cal- 
culations show very little suppression of the quasiparticle 
spectral weight or density of states between the two split 
peaks. A look at Figs. 7c, 8c, and 9e reveals that for 
< 4, N{uj — 0)t w 0.15 and depends very little on 
K^. This is is sharp contrast to the case of antiferromag- 
netic fluctuations discussed in the previous section. This 
difference is to be expected of course, since the antifer- 
romagnetic state is gapped while the ferromagnetic state 
is not. It is clear that none of the diagrammatic approx- 
imations considered here reproduce this spin-splitting of 
the broad quasiparticle peak in ^(k, w) and tunneling 
density of states N{uj) well. In fact, the first order per- 
turbation theoretic result shows a strong suppression of 
the tunelling density of states, which clearly doesn't de- 
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scribe the precursor to the ferromagnetic state well, and 
therefore can be considered to fail qualitatively. We ob- 
serve that the first order perturbation theoretic calcula- 
tion failed to show a suppression of the tunneling density 
of states in the antiferromagnetic case where it is ob- 
tained in the nonperturbative calculations as expected 
(see previous subsection) but does show such a pseudo- 
gap in the ferromagnetic case where it isn't expected and 
doesn't appear in the nonperturbative calculations. It 
therefore qualitatively fails in both cases. Another clear 
sign that not all is well with the perturbation expansion 
is the large quantitative differences between the one-loop 
and two- loop results in Fig. 9, something that could be 
expected at Xz ~ 2.3. In the view of the differences 
between the imaginary time Green's function C/(k, t) ob- 
tained from the Monte Carlo simulations and those of the 
various perturbation-theoretic approximations shown in 
Figs. 9a and 9c which are of the order of O.lt, one would 
expect the spectral functions that produce these rather 
different imaginary time Green's functions to be rather 
different themselves. 



C. Charge fluctuations 

The results of our calculations of the quasiparticle 
Green's function tj(k, r), spectral function ^(k, w) and 
tunneling density of states N{u}) for several values of 
are shown in Figs. 10, 11 and 12. 

For the model studied here, the mass renormalization 
parameter \z is the same for charge and antiferromag- 
netic fluctuations. Therefore the results of the calcu- 
lations for = 24 shown in Fig. 10 correspond to 
\z » 0.05, namely the coupling to the charge fluctua- 
tions is weak. The agreement between the Monte Carlo 
results and those of the various diagrammatic approxi- 
mations is good. As seen in Fig. 10a, the difference be- 
tween the nonperturbative imaginary time Green's func- 
tion and its perturbative approximations at k = (vr, 0) is 
less than the width of the line and of the order of O.OOlt 
for all imaginary times r. If one compares the results of 
the perturbation-theoretic calculations at first and sec- 
ond order, one sees from Fig. 10b that there is virtu- 
ally no difference between the spectral functions A(k, lo) 
obtained by straightforward perturbation theory or the 
self-consistent calculation at either first or second order. 
Hence the slight difference bewteen the first and second 
order calculations, to the extent they aren't an artifact 
of the analytic continuation, must come from vertex cor- 
rections. In contrast to the case of coupling to antifer- 
romagnetic fiuctuations, the spectral functions at second 
order are slightly narrower than their first order counter- 
part. This suggests the first order vertex correction acts 
to reduce the effective charge fluctuation interaction, in 
agreement with the arguments presented in Refi. More- 
over, the nonperturbative A(k, cj) at k = (tt, 0) is slightly 
broader than the second order results, which would indi- 
cate that the higher order diagrams lead to an enhance- 



ment of the effective charge fluctuation interaction, as in 
the case of a coupling to antiferromagnetic fluctuations. 
While this observation is made on the basis of analyti- 
cally continued results, it is consistent with the results 
for other values of presented below, where the en- 
hancement of the effective charge fluctuation mediated 
interaction by higher than second order diagrams can be 
shown to occur on general grounds. 

Fig. 11 shows the the quasiparticle Green's func- 
tion CJ(k, t), spectral function A{k,uj) at k = (7r,0) 
and tunneling density of states N{Ld) for — 4, for 
which Xz ~ 0.45. For this value of corresponding to 
an intermediate coupling regime, the reader will notice 
that the results of the second order perturbation theory 
(self-consistent or not) are not displayed in the figures. 
The reason is that both second order approximations, 
I]^P*(p,ia;„) and E^^^(p,iw„), for the model parameters 
considered here, effectively violate causality requirements 
in that the Eliashberg renormalization factor Z(p,zw„) 
becomes less than one. In terms of the quasiparticle 
self-energy S](p,iu;„), Z(p,«w„) 1 - -^JmT.(p,iuJn) = 
1 _ r+°° du^ irnY.a(pu) ImY.i^i'p.Z) is the imagi- 

nary part of the retarded self-energy and we have made 
use of the spectral representation for the self-energy 

S(p,zc.„) = - '-^ '-"tf^f] . Causality demands 
that the retarded Green's function be analytic in the 
upper-half complex frequency plane and therefore that 
the imaginary part of the retarded self-energy be always 
less than or equal to zero (/mSi{(p, w) < 0) for all val- 
ues of p,ti^. This in turn means that Z(p,ia;„) > 1 
for all values of p,^. One can write the second or- 
der Eliashberg renormalization factor Z'^' (p, iun) = 1-1- 
AZ(i) (p, ia;„)-h AZ(2) (p, ia;„), where AZ(*) (p, iiOn) is the 
change in Z{p,iujn) coming from the i*'' order diagrams. 
AZ*^^^ (p, zwn) is always greater than zero and therefore 
poses no problem as far as the condition Z{p, iLUn) > 1 
is concerned. In the charge-fluctuation case, as was ex- 
plained in Rcf.^, the first order vertex correction has the 
opposite sign compared to the spin-fiuctuation case, and 
leads to a suppression of the effective quasiparticle in- 
teraction. The enhancement of the quasiparticle spin- 
fluctuation vertex comes from the transverse magnetic 
fluctuations that manage to overcome the reduction of 
the effective coupling due to the longitudinal fiuctua- 
tions. Because of this cancellation effect, not only is the 
sign of the first order vertex correction different in the 
magnetic case, it is also smaller in magnitude than in 
the charge-fiuctuation case, under otherwise similar con- 
ditions, as can be seen from the factor 1/3 in Eqs. (|27I28|I 
not present in the corresponding charge-fluctuation case 
in Eqs. H29I30|) . The different sign of the vertex correc- 
tions in the charge and magnetic cases means that while 
in the magnetic case AZ^-^'>{p, iuJn) > and at second or- 
der Z'^^-'fp, zcjn) is always > 1, in the charge fluctuation 
case AZ^^^(p, iw„) < 0. For < 4, we find that the sec- 
ond order contribution to the Eliashberg renormalization 
factor is greater in magnitude than the first order con- 
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tribution, | AZ(2)(p, it^„)| > |AZ(i)(p, iw„)| and since it 
has the opposite sign, Z^-^'>{p,iujn) = 1 + AZ^^'' {p, iujn) + 
AZ^'^\p,iLiJn) < 1- Note that the nonperturbative cal- 
culations always satisfy Z{p,iLUn) > 1, and the problem 
only arises in the perturbative approximation and is a 
sign that, for < 4, the perturbation expansion for the 
charge-fluctuation case is quite badly behaved, possibly 
even more so than for magnetic fluctuations. Aso, the 
fact that the nonperturbative calculations always satisfy 
Z{p, iujn) > 1 is a proof that the higher than second order 
diagrams contribute to an enhancement of the charge- 
fluctuation interaction for these values of k^. 

Fig. lib shows that the spectral function y4(k, w) ob- 
tained from the nonperturbative calculations is notice- 
ably sharper than those produced by the first order self- 
consistent calculations. This means that for g^xo/i = 2 
and Tq — 0.67i the first and higher order vertex correc- 
tions suppress the effective quasiparticle interaction. It is 
therefore not surprising that there are no qualitative dif- 
ferences between the nonperturbative and diagrammatic 
calculations. There are quantitative differences, how- 
ever, and these are more pronounced than in the case 
of a coupling of quasiparticles to antiferromagnetic spin- 
fluctuations for the same value of Xz shown in Fig. 4. 

Finally, Fig. 12 shows our results for — 0.25, for 
which Xz ~ 1.8, hence in the strong coupling regime. 
The difference between the nonperturbative imaginary 
time Green's function and its first order diagrammatic 
approximations seen in Figs. 12b and 12d is a clear in- 
dication of the breakdown of perturbation theory. But 
even in this strong coupling regime, a CDW-precursor 
pseudogap in the spectral function A(k, lu), which can be 
expected to occur on general grounds^ is not seen. The 
pseudogap effects in the charge fiuctuation case thus re- 
quire a stronger coupling still (larger coupling constant 
g'^ or smaller value of k^). 



IV. DISCUSSION 

In Refpi, we showed that the magnetic pseudo- 
gap induced by a coupling to antiferromagnetic spin- 
fluctuations and the spin-splitting of the quasiparticle 
peak induced by a coupling to ferromagnetic spin fluctua- 
tions were not captured by the first order self-consistent, 
or Eliashberg, approximation. The main result of this 
paper, is that these phenomena also lie beyond the two 
magnetic-fiuctuation exchange theories (self-consistent or 
not), which contain first order vertex corrections. While 
this does obviously not constitute a proof, these results 
are consistent with the conjecture expressed in Refi that 
the pseudogap effects found in the Monte Carlo calcula- 
tions are intrinsically nonperturbative in nature. Since 
the calculations reported here show that the first or- 
der vertex corrections alone do not produce a magnetic 
pseudogap, the physics of that state must then mainly 
come from the higher order spin-fluctuation exchange 
processes. The results presented here and in RefJ. also 



indicate that a CDW pseudogap induced by coupling to 
the scalar dynamical molecular field (Eq. must also 
originate from high order charge-fluctuation exchange 
processes. Close enough to a second order CDW tran- 
sition, the diverging CDW correlation length is bound 
to exceed the characteristic length scale for quasiparti- 
cles and the calculations of Ref.^^ showed that when this 
happens a pseudogap opens in the quasiparticle spec- 
trum. The first order vertex correction can't produce 
the pseudogap state, since as we have seen, in the case of 
charge fluctuations it leads to a suppression of the inter- 
action. In fact we even found that for the range of model 
parameters considered here, the second order diagrams 
more than cancel the contribution from the first order 
terms leading to a second order Eliashberg renormal- 
ization parameter Z(p,iw„) < 1, which is inconsistent 
with causality requirements. Moreover, we expect this 
"over-cancellation" effect to get worse as gets smaller 
than the lowest value considered here, = 0.25. Since 
Z(p, iujn) must be > 1 when all the diagrams are summed 
up, as in the Monte Carlo simulations, one can conclude 
that the higher than second order terms must give a con- 
tribution AZ to Z which is positive. Therefore, higher 
order charge-fluctuation exchange processes produce an 
enhancement of the effective quasiparticle interaction, as 
in the magnetic case, and it must be through this en- 
hancement of the effective interaction that a pseudogap 
can appear in the quasiparticle spectrum on the border 
of long-range CDW order. 

These observations lead one to a unified picture of the 
pseudogap state found in our model of quasiparticles cou- 
pled to spin or charge fluctuations. When the dynamical 
molecular field correlation length exceeds the character- 
istic length scale for quasiparticles, either the thermal de 
Broglie wavelengthi2iiii4 or mean free path-, the quasi- 
particles effectively see long-range order and this marks 
the onset of the pseudogap state. This state must be pro- 
duced by high order spin or charge-fluctuation exchanges 
which contain subtle quantum mechanical coherence ef- 
fects. In the magnetic fluctuation case, the first order 
vertex correction favors the pseudogap state, while in 
the charge fluctuation case it suppresses it. This im- 
plies one has to be closer to the boundary of long-range 
charge order to observe a pseudogap than one has to 
be to the boundary of magnetic long-range order, under 
otherwise similar conditions. As the dynamical molec- 
ular field correlation length increases, the mass renor- 
malization parameter Xz gets larger, and therefore the 
many-body effects become stronger. Our results show 
that the agreement between the results of the Monte 
Carlo simulations and the perturbation-theoretic results 
gets worse as Xz increases, and that not surprisingly, 
the perturbation-theoretic calculations break down when 
one enters the strong coupling regime Xz > 1, where 
the pseudogap is found. A more rigourous analysis of 
the relevance of the effective quasiparticle interactions as 
increases or as the energy scales are decreased would 
require a renormalization group (RG) treatment^^. Re- 
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cent RG calculationgi*' on the border of the ferromagnetic 
state indicate that the quasiparticle interactions are in- 
deed relevant in d < 3, and the RG flows to strong cou- 
pling as the energy cutoff is decreased. 

One would like to understand what property of 
the full vertex function ^(x, r; x', r'|x" , t") =< 
Tripa{x,T)ilj\{K' ,t')M'{x" ,t")} > IS responsible for the 
appearance of the pseudogap and seems to be missing in 
the first order approximation to r^ ,^(x, r; x', r'|x" , r"). 

Our physical picture of the pseudogap state emerging 
from quantum mechanical coherence effects contained in 
high order Feynman diagrams is to be contrasted with the 
results of Refsii^iiSiiS where a suppression of the quasi- 
particle tunneling density of states at the Fermi level is 
obtained in the single spin or charge- fluctuation exchange 
approximation. This effect is typically obtained with rel- 
atively large magnetic or charge correlation lengths. In 
our model, the calculations reported here and in Refi 
show that as one approaches the border of magnetic long- 
range order, — > 0, the multiple spin-fluctuation ex- 
change processes become important long before a sup- 
pression of the quasiparticle tunneling density of states 
at the Fermi level is seen in the first order perturbation- 
theoretic and self-consistent calculations. Indeed, pseu- 
dogap effects are only obtined in our calculations when 
the dimensionless mass renormalization parameter Xz > 
1, i.e in the strong coupling regime where one doesn't 
expect diagrammatic perturbation theory to give reli- 
able approximations. The above finding is likely to be 
valid more generally, since the intuitive arguments for 
the physical origin of the pseudogapiiSiiii^ lead one to 
expect the breakdown of Migdal's theorem to be a generic 
occurence near a spin or charge instability. There is also 
an important difference between a vertex correction in- 
duced pseudogap and a single-fluctuation exchange pseu- 
dogap. In the latter case, there is no essential distinc- 
tion bewteen spin and charge fluctuations, in that at the 
single- fluctuation exchange level, for a given fluctuation 
spectrum the spin and charge-fluctuation theories of the 
quasiparticle spectral function can be made identical by 
an appropriate scaling of the coupling constant to the 
molecular fleld. This is no longer the case when vertex 
corrections are included, since these actually depend on 
the nature of the Hubbard-Stratonovich field, in our case 
vector versus scalar. The distinction could turn out to be 
essential, since we find, for a range of model parameters, 
that a pseudogap is observed for quasiparticles coupled 
to spin fluctuations but not in the corresponding charge- 
fluctuation case. 

Moukouri et al^ have developed a many-body theory 
of the precursor pseudogap to the Mott transition in the 
half-filled Hubbard model. Their theory is inspired by the 
fluctuation exchange approximation (FLEX)^° in which 
bare spin and charge susceptibilities are used to build 
up the effective quasiparticle interaction, corresponding 
to g'^xi^T^) ill om' model. The key respect in which the 
theory of Moukouri et aL'^^ differs from FLEX is that the 
coupling to spin and charge fluctuations are not given by 



the bare on-site Coulomb repulsion, but by renormalized 
parameters determined self-consistently in such a way 
that an exact relationship between the single and two- 
particle Green's functions is satisfled. This last step goes 
beyond perturbation theory and it is therefore plausible 
that the precursor pseudogap to the Mott transition seen 
in the Monte Carlo simulations of the half-flUed Hubbard 
model^^ is also nonperturbative in origin. The analog of 
their scheme for the present model would be the use of 
the first order perturbation theory approximation for the 
quasiparticle self-energy S(^^*)(p,a;), Eq. 1^5)1 and a si- 
multaneous renormalization of the coupling constant g 
and the correlation wavevector k^. A renormalization 
of the coupling constant g could account for all vertex 
corrections provided they are local in space and time. 

One can indeed get a pseudogap in the tunneling den- 
sity of states with the first order perturbation theory ap- 
proximation to I](p,ijj) (Eq. (123)), as in Refsiiii^, pro- 
vided is renormalized to lower values and g renormal- 
ized to higher ones. One would thus have to renormal- 
ize the model to stronger coupling, roughly to values of 
Xz ''^ 5. It should be clear that in this regime, first order 
perturbation theory is not controlled. One would also 
naively expect a sensible renormalization scheme that 
goes beyond the one-loop level to lead to renormalized 
values of larger than the bare value. The renormalized 
theory should be further away from the magnetic insta- 
bility than the one-loop approximation rather than closer 
to it, since ideally one would like the improved theory to 
satisfy the Mermin- Wagner theorem in two dimensions. 
If were to be increased by the renormalization scheme, 
in order to obtain a pseudogap in N(lo) one would likely 
need a large renormalization of the coupling g and such 
a scheme for the present model does not look promising 
to the author. However, it is important to note that the 
model studied here, although similar in some respects, 
is actually different than the one considered in Ref.^<^. 
The renormalization scheme proposed by Moukouri et 
al. which works well for the Hubbard model need not 
necessarily apply to other theories. 

In Refi, we pointed out that in the case of quasipar- 
ticles coupled to ferromagnetic spin ffuctuations, our re- 
sults are at variance with expectations based on the stan- 
dard theory of quantum critical phenomena—. Since the 
dynamical exponent z = 3, in d = 2 spatial dimensions, 
the effective dimension is d+z = 5 and is greater than the 
upper critical dimension dc = i above which one would 
expect the first order theory to be at least qualitatively 
correct. But our nonperturbative results show that at 
least for small enough k^, the first order theory quali- 
tatively breaks down. Pepin et. ali& have shown that 
the quasiparticle interactions are indeed relevant in the 
RG sense for ferromagnetic fiuctuations in d < 3, a re- 
sult consistent with our findings. For antiferromagnetic 
and charge fluctuations, d+z — 4, the marginal case, and 
hence the qualitative breakdown of the first order approx- 
imation isn't necessarily inconsistent with the standard 
theory. However, the scaling relations derived in Refiii 
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rely on the applicability of perturbation theory. If the 
pseudogap effects are indeed intrinsically nonperturba- 
tive in nature, a conjecture that is consistent with the 
present work, it opens the possibility that the physics 
in the proximity of a quantum critical point is dom- 
inated by nonperturbative quantum mechanical effects 
and therefore even richer than anticipated in the earlier 
work^^. A number of new ideas in this field have recently 
been proposed^LS^ and a discussion of some fundamental 
problems associated with quantum critical points can be 
found in Rcf.^"^-. 



V. OUTLOOK 

We studied a nonperturbative formulation of the mag- 
netic interaction model, in which quasiparticles are cou- 
pled to a Gaussian distributed dynamical molecular ex- 
change field. Far from the magnetic boundary, the multi- 
ple magnetic fluctuation exchange processes do not bring 
about qualitative changes to the quasiparticle spectrum. 
But as one gets closer to the border of long-range mag- 
netic order, we find, for a range of model parameters, 
that Migdal's theorem doesn't apply and the quasipar- 
ticle spectrum is qualitatively different from its Eliash- 
berg approximation. Moreover, we find that going one 
step beyond the single spin-fluctuation exchange approx- 
imation and including first order vertex corrections, self- 
consistently or not, doesn't help to reproduce the quali- 
tative changes seen in the nonperturbative calculations. 
Near the magnetic boundary, the simple perturbation ex- 
pansion shows signs it is not well behaved, since the sec- 
ond order results differ greatly from their first order coun- 
terparts. The self-consistent, or renormalized perturba- 
tion expansion, which effectively consists in a reordering 
of the diagrammatic perturbation theory, is better be- 
haved in that the differences between first and second 
order are much less pronounced. However, even if the 
renormalized perturbation expansion converges, our re- 



sults show that it is quite likely to converge to the wrong 
answer, which could be explained if the original pertur- 
bation expansion is divergent. 

The intuitive argument for the onset of pseudogap 
behavioriiSiiiii, namely that if the distance quasiparti- 
cles can travel during their lifetime becomes shorter than 
the molecular field correlation length, these quasiparti- 
cles effectively see long-range order, does not explain the 
failure of the single spin-fluctuation exchange approxi- 
mation. As we pointed out in Ref;"'^, one can get in the 
regime where the mean- free path gets shorter than 1/k 
in the Eliashberg approximation, but fail to observe a 
pseudogap in this regime. Since we have not been able 
to produce a good fit to the Monte Carlo simulations 
by including either first order or vertex corrections that 
are local in space and time, i.e by a renormalization of 
the coupling constant g to the molecular field, we con- 
jecture that the physical origin of the pseudogap state 
found in the present calculations lies in non-local ver- 
tex corrections produced by high order spin-fluctuation 
exchanges. These vertex corrections effectively induce a 
quasiparticle coupling to the dynamical molecular field 
that is non-local in both space and time. The above con- 
jecture raise the question of what essential property of 
the vertex function is not captured by its first order ap- 
proximation. A study of the vertex function along the 
same lines as the work reported here for the single parti- 
cle Green's function should provide further insights into 
this problem. 
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J ^ ^ ^ -|- FIG. 2: The various perturbation-theoretic approximations 



^j„,„„^ to the quasiparticle self-energy considered in this paper, (a) 
'\ ,iiii>,. First order perturbation theory, (b) First order self-consistent 
/ \ ^ '\ (Eliashberg). (c) Second order perturbation theory (d) Sec- 
"t" — ^i^>^i^t^ — ond order self-consistent. In this figure a single line de 



notes the bare quasiparticle propagator Go while a double 
line denotes the dressed propagator G to be determined self- 



«•"•», consistently. 
I > i ► I > i ► + 



FIG. 1: Diagrammatic expansion for the single particle 
Green's function. Eg. 12011 up to two spin or charge fluctu- 
ation exchanges. The dashed line represents the dynamical 
susceptibility x(ci, j^^n)- 
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FIG. 3: The various perturbation-theoretic approximations 
to the quasiparticle imaginary time Green's function C/(k, r), 
spectral function A(k, u) and tunnehng density of states N{uj) 
are compared to the results of the Monte Carlo simulations 
for — 24 in the case of quasiparticles coupled to antifer- 
romagnetic spin-fluctuations. 1"* order pt corresponds to the 
approximation to the self-energy shown in Fig. 2a and given 
by Eg. 1251 . 2"'' order pt corresponds to the approximation 
to the self-energy shown in Fig. 2c and given by Eq.l|27|l. 1"* 
order sc corresponds to the approximation to the self-energy 
shown in Fig. 2b and given by Eq.^. 2""^ order sc corre- 
sponds to the approximation to the self-energy shown in Fig. 
2d and given by Eg. 1281 . The error bars on the Monte Carlo 
imaginary time Green's function are not shown for clarity. 
They are of the order of 0.00002t. (a) t/(k,w) at k = (tt.O). 
(b) A{'k,uj) at k= (tt.O). (c) N{lo). 
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FIG. 4: The various perturbation-theoretic approximations 
to the quasiparticle imaginary time Green's function C/(k, r), 
spectral function A(k, u) and tunnehng density of states N{uj) 
are compared to the results of the Monte Carlo simulations 
for K'^ = 4 in the case of quasiparticles coupled to antiferro- 
magnetic spin-fluctuations. 1"* order pt corresponds to the 
approximation to the self-energy shown in Fig. 2a and given 
by Eg. 1251 . 2"'' order pt corresponds to the approximation 
to the self-energy shown in Fig. 2c and given by Eq.l|27|l. 1"* 
order sc corresponds to the approximation to the self-energy 
shown in Fig. 2b and given by Eq.^. 2""^ order sc corre- 
sponds to the approximation to the self-energy shown in Fig. 
2d and given by Eg. 1281 . The error bars on the Monte Carlo 
imaginary time Green's function are not shown for clarity. 
They are of the order of 0.00015t. (a) C?(k,w) at k = (tt.O). 
(b) A{'k,uj) at k= (tt.O). (c) N{lo). 
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FIG. 5: The various perturbation-theoretic approximations 
to the quasiparticle imaginary time Green's function C/(k, r), 
spectral function A(k, u) and tunneling density of states N{lu) 
are compared to the results of the Monte Carlo simulations 
for = 1 in the case of quasiparticles coupled to antiferro- 
magnetic spin-fluctuations. 1"* order pt corresponds to the 
approximation to the self-energy shown in Fig. 2a and given 
by Eg. 1251 . 2"'' order pt corresponds to the approximation 
to the self-energy shown in Fig. 2c and given by Eg. 12711 . 1°* 
order sc corresponds to the approximation to the self-energy 
shown in Fig. 2b and given by Eq.^. 2""^ order sc corre- 
sponds to the approximation to the self-energy shown in Fig. 
2d and given by Eq.JlHJ. The error bars on the Monte Carlo 
imaginary time Green's function are not shown for clarity. 
They are of the order of O.OOlt. (a) g{k,ij) at k = (tt.O). 
(b) A{'k,uj) at k = (tt.O). (c) a(k,tj) at k = (7r/4,7r/4). (d) 
A{k,Lu) at k= (7r/4,7r/4)- (e) N{lo). 
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FIG. 6: The various perturbation-theoretic approximations 
to the quasiparticle imaginary time Green's function C/(k, r), 
spectral function A(k, u) and tunnehng density of states N{uj) 
are compared to the results of the Monte Carlo simulations for 
= 24 in the case of quasiparticles coupled to ferromagnetic 
spin-fluctuations. 1^' order pt corresponds to the approxima- 
tion to the self-energy shown in Fig. 2a and given by Eg. 1)25^ . 
2"'' order pt corresponds to the approximation to the self- 
energy shown in Fig. 2c and given by Eg. 12711 . 1"' order sc 
corresponds to the approximation to the self-energy shown in 
Fig. 2b and given by Eg. 12611 . 2"'' order sc corresponds to the 
approximation to the self-energy shown in Fig. 2d and given 
by Eq.lgHJ. The error bars on the Monte Carlo imaginary 
time Green's function are not shown for clarity. They are of 
the order of 0.00002t. (a) C/(k,tj) at k = (7r,0). (b) A(k,u;) 
at k= (7r,0). (c) N{u). 
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FIG. 7: The various perturbation-theoretic approximations 
to the quasiparticle imaginary time Green's function C/(k, r), 
spectral function A(k, u) and tunnehng density of states N{uj) 
are compared to the results of the Monte Carlo simulations for 
= 4 in the case of quasiparticles coupled to ferromagnetic 
spin-fluctuations. 1^' order pt corresponds to the approxima- 
tion to the self-energy shown in Fig. 2a and given by Eg. 1)25^ . 
2"'' order pt corresponds to the approximation to the self- 
energy shown in Fig. 2c and given by Eg. 12711 . 1"' order sc 
corresponds to the approximation to the self-energy shown in 
Fig. 2b and given by Eg. 12611 . 2"'' order sc corresponds to the 
approximation to the self-energy shown in Fig. 2d and given 
by Eq.lgHJ. The error bars on the Monte Carlo imaginary 
time Green's function are not shown for clarity. They are of 
the order of 0.0002t. (a) C?(k,a;) at k = (7r,0). (b) A(k,a;) at 
k= (7r,0). (c) N{lo). 
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FIG. 8: The various perturbation-theoretic approximations 
to the quasiparticle imaginary time Green's function C/(k, r), 
spectral function A(k, u) and tunnehng density of states N{uj) 
are compared to the results of the Monte Carlo simulations for 
= 1 in the case of quasiparticles coupled to ferromagnetic 
spin-fluctuations. 1^' order pt corresponds to the approxima- 
tion to the self-energy shown in Fig. 2a and given by Eg. 1)25^ . 
2"'' order pt corresponds to the approximation to the self- 
energy shown in Fig. 2c and given by Eg. 12711 . 1"' order sc 
corresponds to the approximation to the self-energy shown in 
Fig. 2b and given by Eg. 12611 . 2"'' order sc corresponds to the 
approximation to the self-energy shown in Fig. 2d and given 
by Eq.lgHJ. The error bars on the Monte Carlo imaginary 
time Green's function are not shown for clarity. They are of 
the order of O.OOOSt. (a) C?(k,a;) at k = (7r,0). (b) A(k,a;) at 
k= (7r,0). (c) N{lo). 
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FIG. 9: The various perturbation-theoretic approximations 
to the quasiparticle imaginary time Green's function C/(k, r), 
spectral function A(k, u) and tunneling density of states N{lu) 
are compared to the results of the Monte Carlo simulations 
for = 0.25 in the case of quasiparticles coupled to ferro- 
magnetic spin-fluctuations. 1"* order pt corresponds to the 
approximation to the self-energy shown in Fig. 2a and given 
by Eg. 1251 . 2"'' order pt corresponds to the approximation 
to the self-energy shown in Fig. 2c and given by Eg. 12711 . 1°* 
order sc corresponds to the approximation to the self-energy 
shown in Fig. 2b and given by Eq.^. 2""^ order sc corre- 
sponds to the approximation to the self-energy shown in Fig. 
2d and given by Eq.JlHJ. The error bars on the Monte Carlo 
imaginary time Green's function are not shown for clarity. 
They are of the order of O.OOlt. (a) g{k,ij) at k = (tt.O). 
(b) A{'k,uj) at k = (tt.O). (c) a(k,tj) at k = (7r/4,7r/4). (d) 
A{k,Lu) at k= (7r/4,7r/4)- (e) N{lo). 
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Charge fluctuations: K = 24 
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FIG. 10: The various perturbation-theoretic approximations 
to the quasiparticle imaginary time Green's function C/(k, r), 
the spectral function yl(k,tj) and tunnehng density of states 
N{uj) are compared to the results of the Monte Carlo sim- 
ulations for K.^ = 24 in the case of quasiparticles coupled 
to charge fluctuations. 1^' order pt corresponds to the ap- 
proximation to the self-energy shown in Fig. 2a and given by 
Eg. 1251 . 1** order sc corresponds to the approximation to the 
self-energy shown in Fig. 2b and given by Eq.lim. The error 
bars on the Monte Carlo imaginary time Green's function are 
not shown for clarity. They are of the order of 0.00002t. (a) 
C;(k,c^) at k= (7r,0). (b) A{'k,uj) at k = (tt.O). (c) N{uj). 
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FIG. 11: The various perturbation-theoretic approximations 
to the quasiparticle imaginary time Green's function C7(k, r), 
the spectral function A(k,Lj) and tunnehng density of states 
N{uj) are compared to the results of the Monte Carlo simula- 
tions for — Am the case of quasiparticles coupled to charge 
fluctuations. 1^* order pt corresponds to the approximation 
to the self-energy shown in Fig. 2a and given by Ea. l|25|l . 1"* 
order sc corresponds to the approximation to the self-energy 
shown in Fig. 2b and given byEq.lgSJ. The error bars on the 
Monte Carlo imaginary time Green's function are not shown 
for clarity. They are of the order of 0. 00021. (a) C/(k,(j;) at 
k=(7r,0). (b) ^(k,w) at k = (7r,0). (c) N{lu). 
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Charge fluctuations: K = 0.25 ; k = (7i/4,7i/4) 



S 





1 1 1 1 

(d) 




1 order pt 




1^' order sc 




Monte Carlo 


... .7. J 


/ > 

/ ' .. 
/ > ■■ ■■. 
/ > .• •• 

■ 1-. . . . 1 . . . 'I 1 . .'.•■.....]-._._■-.■. 


2 


10 12 3 4 



co/t 



FIG. 12: The various perturbation-theoretic approximations 
to the quasiparticle imaginary time Green's function C/(k, r), 
the spectral function A(k,Lj) are compared to the results of 
the Monte Garlo simulations for = 0.25 in the case of 
quasiparticles coupled to charge fluctuations. 1"* order pt 
corresponds to the approximation to the self-energy shown in 
Fig. 2a and given by Ea. ()25^ . 1"* order sc corresponds to the 
approximation to the self-energy shown in Fig. 2b and given 
by Eq.lgSJ. The error bars on the Monte Carlo imaginary 
time Green's function are not shown for clarity. They are of 
the order of 0.002t. (a) g('k,u) at k = (7r,0). (b) A{]s.,uj) 
at k = (7r,0). (c) e(k,tj) at k = (7r/4,7r/4). (d) A(k,tj) at 
k = (7r/4,7r/4). 



